Phase diagram of anisotropic boson t-J model 
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We have studied by Quantum Monte Carlo simulations the low temperature phase diagram of a 
mixture of isotopic, hard core bosons, described by the t-J z -J± model, with J±=aJ z . Separation 
of superfluid hole-rich and insulating, antiferromagnetically ordered hole-free phases is observed at 
sufficiently low hole density, for any a < 1. A two-component checkerboard supersolid phase is 
not observed. The experimental relevance and possible broader implications of these findings are 
discussed. 
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Impressive scientific and technological advances in 
trapping cold atoms in optical lattices (OL) [H, 0, & Q 
render it now feasible to create in the laboratory remark- 
ably close experimental realizations of model many-body 
systems long thought of as of mostly academic interest. 
An example is the isotopic two-component Bose mixture, 
providing a rich playground for many-body physics due 
to the various phases that it is expected to display, in- 
cluding a number of physically distinct superfluid phases 
[H, 0]. The ongoing experimental investigation of this 
system justifies the theoretical exploration of its phase 
diagram, not only to help in the interpretation of exper- 
imental data, but also for the more general purpose of 
guiding the search for novel phases of strongly correlated 
quantum many-body systems. 

In this work, we model a mixture of two different 
species of hard core bosons (the limit of large on-site 
repulsion is assumed for both like and unlike bosons) of 
equal masses, via the two-dimensional (2D) boson t-J z - 
J± model 0: 



H =-t +b\bj+h. 

(ij) V 
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A square lattice of N=LxL sites is assumed, with peri- 
odic boundary conditions (PBC). Two species (A and B) 
of bosons of equal masses are defined, for which a], b], 
are creation operators, whereas hi = a]aj, rhi = b\bi are 
number operators. The sum in (TT]) runs over all pairs of 
nearest- neighboring (NN) sites. 

The Hamiltonian |T]) is defined in the subspace in 
which no double occupation of sites is possible. The 
parameters of the model, namely t, J z and J± are all 
non-negative; henceforth, we shall take t to be our en- 
ergy scale, and set it equal to one. We set J z = J 
and Ji = aJ, a expressing the anisotropy between the 
"antiferromagnetic" coupling J, represented by the sec- 



ond term in (TT]), and the "ferromagnetic" exchange cou- 
pling Jj_, represented by the last term in |T]). In this 
work, we focus on the parameter region in which both 
a and J are less than 1. The hole density is defined 
as h = 1 - (N A + N B )/N, where N A (N B ) is the num- 
ber of particles of species A (B). All throughout, we 
assume N A — Nb < N/2, i.e., the system has no net 
"magnetization". The isotropic version of ([1]) (i.e., with 
a = 1), can be derived from the two-component isotopic 
Bose Hubbard model in the limit of large on-site repul- 
sion and small hole concentration 0, [H, Q; in an opti- 
cal lattice, anisotropy could arise from additional longer- 
ranged (e.g., dipolar) interactions among particles. At 
exactly half filling (i.e., no holes), (JTJ) can also be cast 
in the spin language; for example, it is isomorphic to a 
spin-1 antiferromagnetic Heisenberg model with uniaxial 
single-ion anisotropy, possibly relevant to some magnetic 
systems [Toj | . 

The fermion counterpart of |T]) with ev=l is known as 



the t-J model [ll|, [l2j, ll3[, and has been the subject of 



a wealth of theoretical work, because of its posited con- 
nection to high-temperature superconductivity (HTS). In 
spite of an enormous effort now spanning almost two 
decades the phase diagram of the t-J model remains 
relatively poorly understood. Basic questions, such as 
the presence of a superconducting ground state, are yet 
largely unanswered, essentially due to the lack of a suffi- 
ciently robust theoretical method for strongly correlated 
fermions. On the other hand, the case of Bose statis- 
tics (TT]) can be studied by Quantum Monte Carlo (QMC) 
simulations, yielding essentially exact numerical results. 
For this reason, some studies have used ([T]) as a starting 
point to investigate physical effects such as stripe forma- 
tion, also believed to be relevant to HTS 

The realization of strongly correlated models such as 
(TT|) in OL is the goal of a current, intense experimental 
and theoretical effort. It seems therefore worthwhile to 
provide quantitative, reliable theoretical information on 
phase diagrams of these models, especially to guide the 
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experimental search for novel phases of matter. 

In this work, we study the low temperature phase di- 
agram of if as a function of the hole density h, for 
a < 1 and different values of J. Specifically, we con- 
sider the physically realistic region J < 1 for which the 
isotropic version of (P) displays no phase separation, as 
shown in Refs. 0, LI- Of particular interest is the effect 
of anisotropy on the phase diagram of ([I]), chiefly with 
regard to phase separation, superfluidity, and the possi- 
ble presence of a two-component checkerboard supersolid 
phase (SPS). We use the Worm Algorithm (WA) in the 
lattice path-integral representation [l5| to compute ther- 
mal expectation values of physical operators; the calcu- 
lation requires a relatively straightforward extension of 
the WA, to allow for the simulation of physical processes 
associated with the exchange term of (UJ. 

In the absence of holes (i.e., h = 0), the Hamilto- 
nian |T]) does not differentiate between Bose and Fermi 
statistics; for, upon performing a straightforward basis 
transformation (possible on a bipartite lattice), the SU(2) 
symmetry of the Bose Hamiltonian can be restored. For 
a < 1 , the system features a transition (at finite temper- 
ature T/v) to an antiferromagnetically ordered state, in 
which a particle of type A(B) is preferentially surrounded 
by particles of type B(A) [16[. At a = 1, it is Tn = 0, but 
the ground state is still ordered [ItJ ■ For a > 1, the sys- 
tem is superfluid (SF) below the Berezinskii-Kosterlitz- 
Thoulcss (BKT) transition temperature Tbkt, but with 
no net flow of matter, as the flow of one component is 
exactly compensated by counterflow of the other [18|. 

The presence of mobile holes is generally expected to 
result in a reduction of the antiferromagnetic order, as 
holes scramble it with their motion, as well as to give 
rise to a SF phase of holes at low temperature; of interest 
is the possible existence of a supersolid phase, in which 
both types of orders may coexist in a single homogeneous 
phase. 

The occurrence of checkerboard or antiferromagnetic 
order can be ascertained by numerical simulations of 
the staggered density order parameter x(/i) = \J S(ir, tt) 
where 5(q) = (/5 q /5_ q ), with 
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and where is the position of the ith lattice site, and 
(...) stands for thermal expectation value. The SF density 
Ps(h) of the fluid of holes is computed using the usual 
"winding number" estimator. 

We discuss the isotropic (a=l) case first. As shown in 
Refs. 0, HI) the system features a homogeneous ground 
state, for any hole concentration, for J < 1.5. The AF 
order that exists in the undoped system at T=0, is sup- 
pressed by an arbitrarily small hole concentration. The 
underlying mechanism is simple: holes break the SU(2) 
symmetry of the undoped Hamiltonian in favor of the 
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FIG. 1: Superfluid transition temperature of the uniform hole 
gas in the model |l} with a = 1 (i.e., isotropic), as a function 
of hole density h. Filled circles show numerical estimates of T c 
(in units of t), determined as explained in the text. Dotted 
line at high h is the theoretical behavior predicted by the 
theory of the dilute Bose gas Here, J — 0.4. 



XY-plane and orient the order parameter. We have ex- 
plicitly verified this conclusion by performing finite-size 
scaling analysis of x(/i). This physical result is analogous 
to that observed for a system of lattice hard core bosons 
with a nearest-neighbor repulsive interaction potential of 
strength V — 2t, for which doping away from half-filling 
destroys checkerboard order [20( ■ This points to a signif- 
icant difference between model |T]) with Fermi and Bose 
statistics. In the case of the fermion t-J model the AF 
is reduced but not completely eliminated by hole doping 
2l| since the SU(2) spin symmetry is preserved. 

The ensuing, uniform hole gas is SF at T=0, and fea- 
tures a BKT transition to a normal fluid (NF) at finite 
temperature. The transition temperature T c can be ob- 
tained using the well known renormalization flow and the 
universal jump of the superfluid density, p s , at T c [22j ]. 
Results are shown in Fig. [I] for the case J— OA. 

A richer phase diagram occurs in the anisotropic (a < 
The ground state energy per hole is defined as 



e(h) = 



E(h) - E(0) 
Nh 



(3) 



where E(h) is the total energy of the system in the pres- 
ence of Nh holes. A minimum of e(h) at a finite hole 
density h cr , signals the separation of the system into two 
phases at h < h cr , one with no holes, and the other with 
hole density h cr . 

Fig. [5] shows e(h) as a function of hole density and for 
different values of the anisotropy parameter a. The value 
of J is 0.4. Our numerical results show a minimum for 
e{h) at a finite hole concentration /i cr (a), for any a < 1; 
that is, the system separates into hole-rich and hole-free 
phases for hole doping below h cr (a) (with h cr (a) — > 
as a — > 1). The hole-free phase features AF order (and 
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FIG. 2: Ground state energy per hole e(/i), denned in (|3]), for 
the Hamiltonian (JT| as a function of the hole density h, for 
different values of the anisotropy parameter a (diamonds). 
Boxes show estimates of e(h) for the isotropic (a=l) case. 
Dotted lines are polynomial fits to the data. Numerical cal- 
culations were carried out on a 12x12 lattice. Here, the tem- 
perature T=0.025 t and J = 0.4. 



obviously no hole-based superfluidity), whereas the hole- 
rich phase is SF, but not antiferromagnetically ordered. 
The occurrence of PS in the anisotropic model can be 
understood based on the "string" picture (25|. In the 
anisotropic model, a hole leaves behind, in its motion, a 
string of bosons of either species, displaced by one lat- 
tice site. Thus, separation of the system in hole-rich and 
hole-free phases becomes energetically advantageous at 
low hole density, as a way to limit the damage caused by 
the holes to the antiferromagnetic order. In the a — > 1 
limit, however, quantum fluctuations associated with the 
J± term of "mend" the damage due the hole mo- 
tion, restoring local order [26}. These considerations are 
clearly independent of quantum statistics, i.e., the ought 
to apply to the fermion models as well. 

The separation of the two phases can be visually ob- 
served by examining configurations generated at low T 
by our algorithm on lattices of sufficiently large size, e.g., 
N = 4096 sites (see Fig. [3j) . We have consistently 
observed this effect for J < 1, and for as little as 1% 
anisotropy. Based on this numerical evidence, we argue 
that PS will occur in the T — > limit, at sufficiently low 
hole concentration, for arbitrarily weak anisotropy. 

At finite temperature, entropy favors mixing of the two 
phases and the occurrence of a homogeneous phase. Be- 
cause both SF and AF transition temperatures are finite, 
an obvious question is whether a homogeneous phase fea- 
turing both types of order, namely a "two-component 
checkerboard supersolid (i.e, a superfluid gas of holes in- 
side a checkerboard quantum antiferromagnet formed by 
the two components), may exist at finite temperature. 

Our simulations have not yielded any evidence of such 
a checkerboard supersolid phase. On raising the tem- 



FIG. 3: Snapshot of an instantaneous (/3-averaged) configu- 
ration of the system of a 64 2 -site lattice. Color coding is as 
follows: sites that are drawn in red (blue) are occupied by a 
particle of type A, with a greater probability the bigger the 
size of the circle; sites that are drawn in black are occupied 
by holes and both types of particles with similar probability, 
the smaller the circle, the more likely for the site to be empty. 
Here, J=0.4, J x =0.3, and /i=0.0586 at T=0.025 t. 



perature, the system evolves into either a non-superfluid 
antiferromagnet, or into a SF with no antiferromagnetic 
order. We find that the region where all transition lines 
come close to each other is extremely hard to study since 
all standard finite scaling techniques fail. In Fig. [4] we 
sketch the simplest phase diagram consistent with our 
data, also the one which we find most plausible; in this 
scenario, NF-SF and NF-AF lines meet the coexistence 
dome at one point, close to its maximum. However, at 
this time we can not exclude other scenarios. For ex- 
ample, the AF line may feature a tri-critical point, and 
the BKT transition may be terminated at the first-order 
AFM line. The possible crossing of the second order AF 
and BKT lines, giving rise to a small region of existence 
of a two-component checkerboard supersolid sitting over 
the phase separation dome has been throughly investi- 
gated, but not observed in any of our simulations. 

There are obviously interesting similarities between the 
phase diagram schematically represented in Fig. [3] and 
the basic phase diagram of HTS. In the case of the model 
investigated here, anisotropy is crucial in stabilizing the 
AF phase, which disappears upon doping in the fully 
isotropic model, as shown above. Anisotropy might also 
play a role in shaping the phase diagram of the fermion 
counterpart of (jTJ) , deemed relevant to HTS, as assuming 
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FIG. 4: Computed phase diagram of the boson t-J z -J± model 
on the square lattice. Horizontal axis shows hole concentra- 
tion, vertical temperature. Figure shows actual Monte Carlo 
data obtained for J± = 0.3, J = 0.4, but the same schematic 
phase diagram is obtained for all values of the model pa- 
rameters considered here. Circles represent the normal-to- 
antiferromagnetic transition, triangles the normal to super- 
fluid, whereas boxes show the boundary of the region in which 
coexistence of superfluid and antiferromagnetic phases is ob- 
served. Dashed lines are only meant as a guide to the eye. 
Statistical errors are comparable to the sizes of the symbols. 



a stronger AF coupling in the z direction may be phys- 
ically justified by considerations of interplane exchange 
coupling (27j . As mentioned above, however, in the case 
of Fermi statistics AF order is not completely destroyed 
by the injection of mobile holes. 

The interesting interplay of phases shown in the phase 
diagram of Fig. Q] suggests that an experimental sys- 
tem described by |T]) may be worth investigating in OL. 
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